Non-Gaussian Path Integration 
in Self-Interacting Scalar Field Theories 
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In self-interacting scalar field theories kinetic expansion is an alternative way of calculating the 
generating functional for Green's functions where the zeroth order non-Gaussian path integral be- 
comes diagonal in a;-space and reduces to the product of an ordinary integral at each point which 
can be evaluated exactly. We discuss how to deal with such functional integrals and propose a new 
perturbative expansion scheme which combines the elements of the kinetic expansion with that of 
usual perturbation theory. It is then shown that, when the cutoff dependent bare parameters in the 
potential are fixed to have a well defined non-Gaussian path integral without the kinetic term, the 
theory becomes trivial in the continuum limit. 



I. INTRODUCTION 

In standard model, the scalar Higgs boson is the only 
fundamental particle that has not been discovered yet. 
Our understanding of the scalar sector is based on rather 
indirect observations. For instance, one loop radiative 
corrections in electroweak theory involve Higgs particle 
and the agreement between the theory and the experi- 
ment can be used to estimate the Higgs mass which is 
expected to be less than 225 Gev [1]. Also the scalar 
field is supposed to break the SU{2) x U(l) gauge sym- 
metry down to a U(l) subgroup by acquiring a non-zero 
vacuum expectation value around 247 Gev (determined 
from the experimental value of the W-boson mass and 
the Fermi coupling constant). 

Apart from these and some other similar inputs, the 
scalar sector has not been tested directly. There is no 
experimental data which may be helpful in understand- 
ing the form and the strength of scalar self interactions. 
Specifically, the scalar potential remains to be deter- 
mined. Assuming that the interactions are renormaliz- 
able and tractable in perturbation theory, the potential 
should at most be a quartic polynomial. Further requir- 
ing spontaneous symmetry breaking, the potential can 
uniquely be fixed to have the form A(<fi 2 — /J 2 ) 2 - 

However, perturbative quantization of scalar fields in 
standard model seems unnatural. Namely, there are 
quantum corrections to Higgs mass which are quadratic 
in cutoff A. Therefore, the mass parameter depends very 
sensitively on the short distance behavior of the theory 
and the value of the bare mass should be very fine tuned 
to maintain Higgs mass within a reasonable range. This 
is the so called naturalness or hierarchy problem. 

Moreover, a perturbative renormalization of the self- 
interacting A(f) 4 theory following the ideas of Wilson in- 
dicates that the theory is free or trivial in the continuum 
limit [2]. Although a non-perturbative investigation of 
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the theory on the lattice agrees with naive one-loop per- 
turbative result for renormalization group flow [3], it is 
also possible to argue that in the large N limit there is a 
non-trivial fixed point in the four dimensional scalar field 
theory [4]. 

Therefore, it is of great importance to develop new an- 
alytical methods alternative to usual perturbation theory 
which may help us to resolve these and similar issues in 
scalar field theories. In this paper, we consider one such 
framework, known as kinetic (or strong coupling) expan- 
sion, where the interaction potential, which can be a non- 
polynomial function, is treated exactly. The main strat- 
egy here is to deal with the kinetic term perturbatively 
and carry out Euclidean path integration in x-space (de- 
fined with a cutoff), where the potential is "diagonal" 
and the functional integral reduces to the product of an 
ordinary integral at each point (see e.g. [6]- [16]). 

In this work we first discuss how to calculate a non- 
Gaussian functional integral which appears in the zeroth 
order contribution in the kinetic expansion. It turns out 
that the bare parameters in the potential should be taken 
cutoff dependent to get a finite result (see e.g. [7]). How- 
ever, the higher order corrections in the expansion be- 
come ill defined due to the presence of uncontrollable 
infinities. One can try to perform an extrapolation to 
cure the problem as proposed in [10, 12], however this 
procedure likely fails in general as discussed in [11, 14]. 
The main obstacle in this program is that there is an 
expansion in derivatives about a zeroth order configu- 
ration which is calculated without paying attention to 
locality since the derivative (kinetic) term in the action 
is removed [15]. 

In section III, an alternative expansion scheme is pro- 
posed where one can overcome this obstacle. We find 
that by introducing a Lagrange multiplier it is possible 
to invert the kinetic term as in a free theory while the po- 
tential can also be integrated out non-perturbatively (but 
order by order in cutoff) yielding a complementary series. 
In section IV, we use this result to show that if the cutoff 
dependences of the bare parameters in the potential are 
chosen specifically to have a well defined non-Gaussian 
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path integral without the kinetic term, then the theory 
becomes trivial in the continuum limit. 



II. PERTURB ATIVE EXPANSION OF THE 
KINETIC TERM 

In this section, we mainly review the kinetic expansion 
of the generating functional (see e.g. [10]). Consider a 
self-interacting, real scalar field theory defined in n di- 
mensional space-time. The Euclidean path integral for 
the generating functional can be written as 



where 



Z[J]=N J D^eU 13 ^- 



v{4>) + J<t>\ 



(1) 



where J means integration over d n x (the x dependence of 
functions are suppressed) N is the (infinite) normaliza- 
tion, P is a constant and V{fa) is the interaction potential 
including a possible mass term. We take n > 1, since we 
will freely play with the constant TV in (1). This is not 
possible when n = 1 since in quantum mechanics the in- 
tegration measure is uniquely fixed by the normalization 
of the state vectors. 

We evaluate (I) by expanding the exponential of the 
kinetic term so that (3 is the perturbation parameter. The 
zcroth order contribution is 



Z [J] = N D<j>e 



I 



\J -V(4>) + J4>] 



(2) 



We first approximate such functional integrals by divid- 
ing n-dimensional space into regions of volume e. One 
can think that the theory is defined with 



e = L n , A = l/L, 



(3) 



where L is the spatial lattice size and A is the cutoff in 
momentum space. In this approximation 
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and path integral measure is given by 
D<p — » ITj dfa, 



(4) 
(5) 

(6) 



where the index i runs over the infinitesimal regions. The 
Dirac delta function 5(x — y) and the derivative operator 
can be represented as 5ij/e and (d(f>)i — (fa+i — fa)/e- 
We will keep e infinitesimal but finite and let e — > to 
reach the continuum limit. 

Eq. (2) can now be rewritten as 

Z [J] = N J Ylidfa e E,[- v W')+J«*i]« 

= NIU Zi(Ji), (7) 



r 

Zi(Ji) = / 



dfa e 



-V(<p t )e+ Ji4.it 



(8) 



Thus the problem reduces to the calculation of an ordi- 
nary integral. Although for simple polynomial potentials 
the exact result of (8) can be found, 1 for our purposes a 
series solution in Ji is enough. Of course the interaction 
potential should ensure the finiteness of each term and 
the convergence of the sum in the series. We take V to 
be an even function V(— <j>) = V{<j>) with V(<f>) — ► +<fr p 
for some p > 1 as (f> — > ±oo to satisfy this technical con- 
dition. Let us scale fa — ► fa/y/e in (8) to organize the 
expansion. It is important to note that any overall mul- 
tiplicative constant (like 1/y/e which appears after the 
scaling) can be absorbed in the definition of the constant 
N in (7) and therefore ignored. Introducing 



and 



V{fa) = V(fa/sfe)e, 



a m = — - / e- y (") u m du, 
ml 



(8) can be evaluated as 



Zi 



(j,) = i + E — J f m - 

m=l 



(9) 
(10) 

(11) 



where ao is also factored out to ensure the normalization 
Zj(0) = 1. Note that only even powers of Ji survive since 
we take V to be an even function. 
Using (II) in (7) one obtains 



z [J] = n, 



= exp 



a a 

fl 2 T 2 



0(e), 



and in the continuum limit we get 



Z [J] = exp 



02 

a 



J{x) 2 d n x 



(12) 



(13) 



Although (13) resembles a functional obtained after a 
Gaussian path integration, there is a crucial difference; 



For V = \<t> 4 the result of (8) is 
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where p F q [a; b; x] is the generalized hypergeometric function 
which is regular at x = and normalized so that p F q [a; b; 0] = 1. 
A standard integral table can be found in [5] 
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there are actually 0(e) corrections to (13) which can be 
seen in (12). These can be ignored at this point, however 
they are not negligible when (3 corrections are calculated 
(see below). 

The constants do and a 2 might depend on the cutoff 
so there can still be infinities hiding in Zq[J]. To get a 
finite result in the continuum, V defined in (9) should 
not depend on e as e — > 0. This can be ensured by taking 
the bare parameters in the potential cutoff dependent. 2 
For instance, if 



V = m 



2 J.2 



+ A< 



(14) 



then for V to become cutoff independent A should be 
fixed as 



(p-2)/2 

A = I - ) A. 



>(p-2)/2 



(15) 



where e M corresponds to a fixed scale n, i.e. = fi~ n , 
and A M is the "renormalized" coupling constant. Using 
(15) the cutoff A and bare coupling A dependences of 
ao and ai can be traded with \i and A M , respectively, 
which remain finite in the continuum limit. Note that for 
the bare mass m there is no need to impose any cutoff 
dependence. Eq. (15) can be viewed as the analog of the 
renormalization group flow in this framework. For p > 2 
the bare coupling vanishes and for p < 2 it diverges in 
the continuum limit e — > or A — > oo. 

To illustrate the above formulas with a specific exam- 
ple, let us take V — m 2 <fi 2 + A</> 4 . Then the constants ao 
and ci2 can be found using (10) which gives 



a 



a.2 = 



(16) 



4 U, 

2 

m 
T 



3/4 



r[3/4] ^[3/4, 1/2; a;] 



A, 



5/4 



T[5/4] !^ [5/4, 3/2; x], (17) 



where x = m 4 e M /(4A AJ ), K is the Bessel and p F q [a;h;x] 
is the generalized hypergeometric function (see also foot- 
note 1). From these, it is clear that A M dependence of 
(13) is not analytic and cannot be obtained in usual per- 
turbation theory. 

The full generating functional (1) can formally be ex- 
pressed as 



Z[J] = exp 







J G(x,y) 



S 2 



5J(x)6J(y) 



d n xd n y 



Z [J] 
(18) 



2 Recall that in usual perturbation theory cutoff dependences of 
the bare parameters in the Lagrangian are fixed by requiring 
finitcness of physical quantities. Here we employ the same phi- 
losophy. 



where G(x,y) = U x S(x — y) and Zo[J] is given in (7). 
Eq. (18) is ready for a perturbative expansion in /3 where 
the zeroth order contribution is given by (13). It is im- 
portant that one uses (7) (together with (11)), but not 
simply (13), for Z [J] in (18) to calculate higher order 
corrections. This is because, from (5), the functional 
derivatives have a (singular) cutoff dependence and when 
acting on 0(e) corrections in (12) this would yield non- 
vanishing contributions in the continuum limit (this is 
actually the case as we will see). 

Expanding (18) to first order in /3 and using (7) one 
can find 



Z[J] = Z [J] + f3 Z [J] J2 G i 



din Zi dhiZn 



dJi dJj 



+ pz [J]Y,Gu^jP + O(p 2 ) (19) 

i 1 

where denote the entries of G(x, y) on the lattice. 
Note that e 2 factor coming from the integration measure 
d n xd n y in (18) is canceled by the cutoff dependence of 
the functional derivatives (5). 

Let us now take the continuum limit. Using Zi given 
in (11) we see that the leading order e terms in (19) 
are e 2 and e for the first and the second lines, respec- 
tively, which are exactly the required powers to convert 
sums into integrals. Unfortunately, however, the diago- 
nal entries of Gij diverge for the kinetic term since in the 
continuum we have G(x,y) — O x S(x — y), so there is a 
possible divergence problem in the second line in (19). 

For the moment, assume that G(x, y) is a smooth and 
finite function. Then in e — > limit (19) becomes 

Z[J] = (1 + /3C) Z [J] 



+4/3 
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Z [J] J G(x,y)J(x)J(y)d n xd n y, (20) 



where C = 2(a 2 /a ) / G(x,x) d n x is a constant. The 
properly normalized generating functional up to 0(/3 2 ) 



is 



Z[0] 

+4/3 



- Z Q [J] 



«2 
«0 



Z [J] G(x,y)J(x)J(y)d n xd n y, (21) 



so the dependence on a (possibly infinite) contribution C 
drops out (not that we set Z a [0] = 1). Of course one can 
now use (13) in (21) since there is no functional derivative 
acting on Z [J]. The constant C is related to "vacuum 
to vacuum" amplitude and it may also be eliminated by 
redefining the normalization N in (7). Thus, the first 
order result (21) is finite in the continuum for smooth 
G(x,y). 

To see what happens for a singular function let us take 
G(x,y) = —m 2 S(x — y). Since the mass term can actu- 
ally be treated exactly, the perturbative first order result 
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may be derived from the exact formula by expanding the 
coefficient of J(x) 2 in (13) around m = for the poten- 
tial V = Va + m 2 <j) 2 . Using, on the other hand, the lattice 
approximation Gij = —m 2 5ij/e in (19) one can directly 
get 



z[o] - 

+2 (3 m 2 



Zo[J] 



604 
a 



Z Q [J] / J{xfd n x. (22) 



As in (20), there is a CZo[J] term in Z\J\ which is elim- 
inated after dividing by Z[0] where the constant C di- 
verges like C ~ J d n x. 

For the kinetic term the situation is more compli- 
cated. Using the approximations for the delta func- 
tion and the derivative operator given below (6) one 
finds Gij = (Si+2,j — 2<Jj + ij + Sij)/e 3 . Therefore 
Gu — 1/e 3 (no sum in i is implied). However, it 
is also possible to define an asymmetric derivative op- 
erator such as d(f>i — (<f>i+i — 0,_2)/(3e) which gives 
Gij = (^ + 2,j-2-2(5i_ij+^_4,j)/(9e 3 ) and thus G li = 0. 

Note that (d<f))i = (d(f>i + d<j>i-\ + d<fii-2)/3 so d in- 
volves an extra averaging over neighboring lattice points 
in terms of d. In the continuum limit both d and d 
asymptote to the same differential operator. But, one 
should actually specify the function space in the contin- 
uum more precisely. Namely, d can only be used as the 
derivative operator if the test functions are not changing 
appreciably along three neighboring lattice points. 

To resolve the issue we use the Fourier transform with 
a momentum cutoff A to specify smoothness of the field 
variables in the continuum. From 



G(x, y) = - J e ^-2V d n p, 



(23) 



one gets G(x,x) <~ — J p n+1 dp <~ — A"+ 2 and thus we 
take Gu = ce -(™+ 2 )/™ in (19) where c is a (finite) numer- 
ical constant. In the continuum limit we then have 



Z[J] 



Z[0] 

-6/?ce- 2/ ™ 



Z [J] + 4/3 



2ct4 

a 



02 



Z [J] J J{x)U x J{x)d n x 
Z [J] J J{xfd n x. (24) 



Unfortunately, the second line of (24) diverges like 
e - 2 /™ ~ A 2 . However, from (22), this divergence can be 
canceled out by a mass counterterm (5m) 2 ~ (3 A 2 . So 
we conclude that, to first order in (3, only a mass renor- 
malization is required after which all Green's functions 
become finite. 

As for the higher order (3 corrections, it is not difficult 
to obtain the analog of (19) on the lattice. For smooth 
and finite G(x,y), 0(/3 2 ) contribution turn out to be fi- 
nite in the continuum. However, analyzing the result for 
the kinetic term to this order we find that some Green's 
functions diverge (e.g. quadratically in A) and it seems 



these cannot simply be removed by local counterterms as 
in (24). We also expect more and more infinities to show 
up in higher orders coming from the multiplications of 
the singular distribution U x 8(x — y). So it seems that 
new regularization or renormalization techniques should 
be developed for higher orders. 3 As pointed out in the 
introduction the main problem here is that in calculating 
Zq[J] locality is neglected since the kinetic term in the 
action is removed. On the contrary, the higher order cor- 
rections are sensitive to local variations since more and 
more powers of the Laplacian appear in the expansion. 



III. AN ALTERNATIVE EXPANSION 

Our main strategy in this section is to use a Lagrange 
multiplier in the path integral to separate the kinetic and 
the potential terms so that (1) can be rewritten as 

Z[J} = J Deb D$ DpeU? <^>- V W + + *(*"*)] 

(25) 

where the normalization N is suppressed. Evidently in- 
tegrating over p gives a delta functional and a further 
trivial integration over $ gives (1). However one can 
now carry out the Gaussian <p integration to get 

Z[J] = f D$Dp eU ^(^-'ri-vW+JW-*)] 

where we have introduced an auxiliary field K. Since 
differentiating with respect to K(x) gives ip(x), the last 
equation can be expressed as 



Z\J\ = cxp 



6 



-d n x 



4/3 J 8K{x) x 5K(x] 
J D<5>DpeU- v ^ +J ' s ' + ^ K -^] 



K=0 



At this point, p and <& integrations become trivial where 
the first one gives a delta functional and the second one 



3 One possibility is to set up a renormalization group flow starting 
from IR. For instance, from (23) one may decompose 



G(x, y) = G(x, y, n) + G(x, y) , 



where 



G(x,y,n) 



J —fl 



ip(x-y) p 2 d n p ^ 



fi is a fixed IR length scale and G contains the integral in (23) 
starting from fi to oo. Now, in (18) one can first calculate the 
contribution of G(x,y, fi). Since G(x,y,fi) is non-singular the 
exact result of the computation can in principle be performed 
to any given order without encountering any infinities or am- 
biguities. This gives a new functional Z^[J] and one can now 
integrate the modes from fi to 2fi to obtain Z2fi [J] etc. 
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sets $ = K. After these integrations one finds 



Z[J] 



5 



-1 _S_ 

SK J 



,[/-V(Jf) + JJf] 



(26) 

which can be evaluated order by order in 1/(3 by expand- 
ing the first exponential. Comparing to (18) where (3 is 
the expansion parameter, we see that (26) gives a com- 
plementary series, reminiscent of a strong- weak coupling 
duality. 

It is not difficult to show that for polynomial potentials 
(26) is equivalent to usual perturbation theory. Noting 
that the functional derivative of exp(J JK) with respect 
to J(x) (K(x)) gives K(x) (J(x)), (26) can be rewritten 
as 



Z[J] = e-^f&^^ef-^ef JK 



— e-f v ^ " - 1 e 4 / 



- 1 _L- 

sif e 



if=0 



= e I- v ^e-^I Jn ~ lj , 



K=0 



(27) 



where in the second line we have used the fact that J and 
K derivatives commute and in the last line exp(J JK) is 
dropped since any J derivative acting on it vanishes after 
setting K = 0. This proves the equivalence since (27) is 
nothing but the usual perturbative expansion formula. 
From the commutator 



/ 



J-O- 1 — K(x) 
SK SK' [ > 



2n: 



6K(x) ' 



(28) 



one can also directly verify that (26) obeys Schwingcr- 
Dyson equation 



2/3 n x 



5Z[J] 
5J{x) 



V 



SJ(x) 



Z[J] + J(x) Z[J] = 0, (29) 



where V = dV/d<j). Note that in (26), for a fixed order in 
1/(3 (and for polynomial potentials), only a finite number 
of terms survive if one expands exp(— J V(K)) since one 
sets K — at the end. 

Eq. (26) is preferable over (27) for organizing the per- 
turbative expansion in powers of 1/(3. Another advantage 
is that non-polynomial potentials can also be treated nat- 
urally. Assuming that V{K) is an infinitely diffcrentiable 
smooth function, (26) shows that there is an interaction 
vertex for the p'th derivative of V{K) if (PV/dRP ^ at 
K = 0. One may reach the same conclusion in perturba- 
tion theory by expanding V(K) around K = 0. 

Having rederived the usual perturbation theory in a 
different way, we now proceed by obtaining another ex- 
pression for the generating functional. Let us start with 
(25) in which the external current J is coupled to <f> rather 
than Integrating over $ one gets 

Z[J] = [ D4>DpZ [ip]eU 'HW)+J* + ip{K-4>. 

J K=0 



where the functional Zq is defined in (2) and we again 
introduced an auxiliary field K. Using the same trick we 
employed below (25), we find 



Z[J] = Z 



SK 



y [J H(KUK)+.Jk] 



(30) 



Adding and subtracting J jn _1 J/(4/3) term in the ex- 
ponential a complete square can be obtained. Shifting K 
as K(x) — > K(x) — U^ 1 J{x)/{2(3) one finally reaches 



z[j] = e \r*el ' Ja ~ lj ] z {) 

or more conveniently 

Z[J]=e[-^ Jn ~ lj \ Zo 



5K 



KUK\ 



2f}UK=J 

(31) 



S ' 
n 6J 



(32) 



where we set (3 = 1 /2 since it can no longer be used as an 
expansion parameter (note 1/(3 factors in (31)). Unlike 
the expression given in (18), we now manage to invert 
the kinetic term as in a Gaussian integral. Using the 
identity (34) below it is easy to verify (32) in the free 
theory when the potential contains only a mass term. 
Note that apart from the first exponential (and a sign in 
the third one) logZ o [D0] can be viewed as an effective 
"potential" added to a free massless scalar field. 

In the continuum limit, the functional Zq is given in 
(13). However, 0(e) corrections are not negligible as be- 
fore. Therefore, a new expansion scheme can be obtained 
by evaluating Z order by order in e and using the result 
in (31) or in (32). 

Let us proceed by a direct computation of the first two 
terms in this expansion (a systematic treatment will be 
presented in the next section). From (12) we obtain 



Zo[f] = exp 

«4 



/ 



1 



K € 



«0 



«2 

2a 2 



f(x) 2 d n x 



,f{xfd n x + 0{n 2 ) 



(33) 



where n is a formal expansion parameter which counts 
the powers of e (more comments on this below). To cal- 
culate the contribution of the exponential wc note the 
identity 



; [ / 2 7J? ] g [/ 2B ] _ „ [/ 2(A + B) ] 



(34) 



which can be verified by setting up a simple Gaussian 
path integral. In (34) , an overall irrelevant (infinite) mul- 
tiplicative constant is omitted on the right hand side. Us- 
ing (33) and (34) in (32) the zeroth order term can be 
found as 



where 



Z\J\ = e 

Pi(x,y) = 



^/ JPlJ ]+0( K ), 



□ 



5{x-y). 



(35) 



(36) 
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The same equations show that to calculate the next order 
contribution one confronts the term 




5 



x 5J(x) 



d n x 



exp 



ju~ l p»j 



where 



— Pi. 

2a 2 



, (37) 



(38) 



The result of (37) contains the x-integral of 

3(aP 2 )L 

+ 6(nP 2 ) X x J P2 (x, yi)P 2 (x, V2)J(yi)J(y2)dy 1 dy 2 

+ J P 2 (x,y 1 )P2(x,y 2 )P 2 (x,y 3 )P 2 (x,y 4 ) x 

J(yi) J(V2) J{Vs)J{yi) dyi dy 2 dy 3 dy 4 . (39) 

The last line can be dropped in the continuum since, from 
(33), the whole expression is multiplied by e. For the first 
two lines, one needs the numerical value of (OP 2 ) xx which 
can be found from the Fourier transform 



□P 2 (x 



■»>~/r 



+ 2p 2 a 2 /a Q 



Mx-y) d n p (4Q) 



When x = y the integral diverges like A™ and thus 
(UP 2 ) XX = c/e where c is a finite constant. The first line 
in (39) yields an infinity, however this can be pushed to 
order k 2 after dividing Z[J] by Z[0]. On the other hand 
1/e factor coming from (OP 2 ) xx in the second line is can- 
celed by the e term in (33), thus there is no problem in 
taking e — > limit. Combining all these, the continuum 
result becomes 



Z[J] 

J P2 (x , yi ) P 2 (x, y 2 ) J(yi ) J(y 2 ) dx dy x dy 2 + O (k 2 ) 

(41) 

which does not contain any infinities. 

Higher order k corrections can be calculated in a sim- 
ilar fashion where one now encounters more coincident 
functional derivatives acting in (32) and there are both 
connected and disconnected contributions. Of course, 
more coincident derivatives give a much more singular 
behavior. However, these are multiplied by higher pow- 
ers of e and there is a chance that one gets a finite result 
in the continuum as for the first order calculation above. 
This is indeed the case as will be discussed in the next 
section. 

The constant k is just a formal expansion parameter 
and one should actually set k = 1 at the end. There- 
fore, the higher order contributions might in principle be 
greater than the lower order ones. However, by analyzing 
the expansion of Zq we observe that higher order terms 



contain coefficients with larger k. Moreover, at least 
for polynomial potentials, a scaling argument shows that 
the coupling constant dependence of a,k has a hierarchy 
fixed by k i.e. for larger k one gets more powers of the 
inverse coupling constant (see e.g. the expressions given 
in (16) and (17)). Therefore, the expansion in n can be 
reorganized as an expansion in the inverse coupling con- 
stant. 



IV. TRIVIALITY 

In section II, we found that to have a well defined 
non-Gaussian path integral in the kinetic expansion the 
bare coupling constant A should obey (15). Therefore, A 
vanishes for p > 2 in the continuum limit as A — > 00, and 
one may wonder whether the theory becomes free in this 
limit. Although the interactions vanish classically, this is 
not an obvious question in quantum theory. Indeed, there 
is no sign of triviality in the kinetic expansion studied in 
section II. 

Let us try to analyze the situation in usual perturba- 
tion theory with the interaction potential (14) (so p is 
assumed to be an integer). Without loss of of generality 
we take p > 2 and focus on the connected graphs. It 
is well known that infinities generically arise in evaluat- 
ing Feynman diagrams, however in our case these can be 
suppressed by the bare coupling constant since it van- 
ishes like (15) as one removes the cutoff A. Consider a 
connected graph with / internal lines, E external lines, N 
vertices and L loops. We have the following topological 
identities: 



21 + E = Np, 



L = I-N + 1. 



(42) 



The superficial degree of divergence D of this diagram 
is given by (for each loop the momentum space integra- 
tion measure contributes n units, and each internal line 
propagator adds —2) 



D = nL- 21. 
Using (42) to eliminate L and I one gets 



(43) 



nE „ T Nn(p-2) 
D = -—-2I + n+ ^ >-. (44) 

The order of the diagram is N and the degree of sup- 
pression S coming from the cutoff dependence of A in 
(15) can be found as 



Nn(p - 2) 



(45) 



Comparing the two numbers we find that S > D . There- 
fore, as A — > 00, \ N term becomes more singular than 
the outcome of the momentum space integral and the 
final result for the diagram should tend to 0. Accord- 
ing to this naive power counting, the only non-vanishing 
graph is the tree level two point function, i.e. the quan- 
tum theory is actually free. It is important to emphasize 
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that this is not a formal proof of triviality since the cal- 
culated degree of divergence D arises only from regions 
of momentum space in which all internal n-momenta go 
to infinity together. It is known that additional diver- 
gences can also show up when the momenta belonging to 
some subgraph tend to infinity. Therefore, D does not 
necessarily give the actual degree of divergence of a dia- 
gram. However, the above power counting is suggestive 
and indicates triviality. 

Let us now try to use the alternative expression (31) 
in evaluating the generating functional. From (12) one 
finds 



z o[f] = cxp 



a 



x Z mt [f], (46) 



where Zi nt [/] has an expansion of the form 



Zint[f] = CXp 



OO r. 

]Te fe -M fe / f(x) 2k d n 

k=2 J 



(47) 



and Ak arc finite coefficients which can be determined in 
terms of dfc, e.g. A 2 = a^ja^ — a\j{2a^). Using (46) and 
(34) in (31), we get 



Z[J] = exp 



J int 



2 

_S_ 

Jk 



UK=J 



where 



A{x,y) = UP 2 (x,y). 



(48) 



(49) 



From (47) the second line of (48) is equivalent to the 
generating functional of a scalar field ip having the inter- 
action potential Vint, 



V ir , 



OO 

k=2 



2k 



(50) 



and the propagator A (here K plays the role of the exter- 
nal current coupled to this hypothetical scalar p) . There- 
fore, it can be evaluated order by order using the well 
known perturbation theory techniques, i.e. by calculat- 
ing the Feynman diagrams corresponding the interaction 
potential (50). The only difference is that an extra 
factor should be attached to each external line since one 
should set UK = J at the end. 

Consider a connected graph in p theory with / internal 
lines, E external lines, Nk vertices of type <p> 2k and L 
loops. For this diagram we have 

21 + E = ^JV fc (2fc), L = I-N + 1, (51) 

k 

where N = ^2 k Nk is the total number of interaction ver- 
tices. Upto a finite constant prefactor, the corresponding 
truncated function becomes 

N I 

G~Snn / d n Pl S(m)A( Pl ), (52) 

m=l 1=1 J 



where 5(m) represents the Dirac delta function at the 
vertex m, A(p) is the momentum space propagator which 
can be calculated from (40) and (49) as 



A(p) = 



p- 2 + 

r o 



(53) 



and S is the contribution of the cutoff dependent coupling 
constants in (50): 



(54) 



Note that the asymptotic behavior of A is very different 
than the propagator of a conventional scalar field and we 
have 



(55) 



Using this inequality in (52) (and ignoring the finite 
parts) one finds 



n i ,. 

g < s n n / d > i ^ m )- 



(56) 



The delta functions impose linear relations among the in- 
ternal momenta which reduce the number of independent 
momentum variables to the number of loops (of course 
one delta function survives which ensures conservation 
of external momentum). Moreover each loop integral in 
(56) diverges like A" and thus 



G < SA 



nL 



(57) 



which gives an upper bound for the actual degree of di- 
vergence of this diagram. Solving for / and L from (51), 
and using (54) one finally obtains 



G < K n( - 1 - E / 2 \ 



(58) 



Eq. (58) shows that all graphs which have more than two 
external lines vanish and all corrections to the propaga- 
tor are bounded in the continuum limit. Therefore, the 
theory becomes trivial as one removes the cutoff. 

A few comments are now in order. Although correc- 
tions to the propagator are all finite, it is not clear if the 
sum of these terms would converge or if the final result 
would correspond to a physical particle. The first term 
in this series was actually calculated in (41). This cor- 
responds to a one-loop diagram in p theory which arises 
from the first term in the potential in (50) (see figure 1). 
The tree level contribution of the same term to the four- 
point function is suppressed by the cutoff dependence of 
the "coupling constant" in Vi n t ■ Note that in the second 
line of (41) the propagator P 2 is attached to the exter- 
nal lines instead of A. This is because, as pointed out 
above, the external lines should be multiplied by in 
obtaining Z[J] and we have n _1 A = P 2 . 
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FIG. 1: The first order correction to the propagator which 
arises from eA2 <p 4 term in the potential (50). The loop in- 
tegral, which is calculated with the propagator A in (53), 
diverges like A n . This divergence is compensated by the e 
factor and the result becomes finite in the A — > oo limit. 

The above argument of triviality also works for non- 
polynomial potentials which cannot be analyzed in usual 
perturbation theory. The only requirement is that the 
cutoff dependences of the bare parameters should be fixed 
as in (15) to have a well defined non-Gaussian path inte- 
gral without the kinetic term. In general (15) is not valid 
for the bare parameters obtained in usual perturbation 
theory; for instance in A0 4 theory in 4-dimensions there 
is a one loop logarithmic A dependence in A. So, the 
triviality argument is not directly applicable to pertur- 
batively treated self-interactions. 

Returning back the kinetic expansion discussed in sec- 
tion II, one may conclude that the infinities arise simply 
because locality is ignored in calculating the zeroth or- 
der contribution Z [J] and, oppositely, the corrections 
demand locality (this viewpoint was emphasized in [15]). 
Therefore, to be able to use (18) one my try to find a 
way of resolving this dilemma. One possibility is to set 
up an RG flow starting from IR as discussed in footnote 
3. It would be interesting to verify triviality directly from 
(18) by using this RG flow, so that the locality can be 
recovered step by step starting from IR. 

V. CONCLUSIONS 

In this paper, an alternative expansion scheme is ob- 
tained in self-interacting scalar field theories which uses 
the elements of kinetic expansion and usual perturbation 
theory together. The central object in this framework 
is the non-Gaussian functional integral of the potential. 
These integrals naturally arise in the kinetic expansion, 
which, however, is plagued by the presence of infinities. 
Kinetic expansion is ill defined since locality is completely 
neglected in the zeroth order configuration around which 
there is a local expansion. To resolve the problem we 
introduced a Lagrange multiplier to separate the kinetic 
and the potential terms in the path integral. In this way, 
we are able to invert the kinetic term while the potential 
can be integrated out non-perturbatively but order by or- 



der in a cutoff parameter yielding a series expansion for 
the generating functional. Although we just consider a 
real scalar field, it is not difficult to generalize the formu- 
las for coupled multiple scalars. For instance, the analog 
of the one dimensional integral (8) is a multi-dimensional 
one fixed by the number of scalars in the theory. 

Using this new scheme, we show that if the bare pa- 
rameters in the potential are chosen to have a well de- 
fined non-Gaussian path integral then all graphs except 
the two-point function vanish in the continuum limit. 
Therefore, the theory becomes trivial when the cutoff 
is removed. Note that the bare couplings are fixed "kinc- 
matically", i.e. to define a proper integration measure. 
For polynomial potentials it should be possible to give 
a proof of triviality in usual perturbation theory as sug- 
gested by the naive power counting discussed in section 
IV. However, it is difficult to see triviality in the kinetic 
expansion. 

Changing the potential only modifies the constants a m 
(i.e. their dependences on the parameters of the theory) 
and the other steps in the calculations remain unaffected. 
This indicates that the self-interaction potential simply 
plays the role of a weight function in the path integral. 
Note that in Euclidean space the convergence of the gen- 
erating functional (1) should be more rapid for "larger" 
interaction potentials since contribution of each field con- 
figuration in the functional integral would decrease due 
to the exp(— J V) factor. This can be observed in our 
approach since the path integral containing the potential 
reduces to product of ordinary integrals. However, this 
is contrary to what we see in usual perturbation theory; 
for example in four-dimensions although A</> 4 theory is 
perturbatively tractable A</> 6 is not. 

It would be interesting to extend the present work in 
different directions. Firstly, it is possible to combine 
the above formalism with the usual perturbation theory. 
This can be achieved by taking all bare quantities in the 
Lagrangian (the wave-function, the coupling constants 
etc.) as provided by the perturbative renormalization 
theory and use them in the expression (31) or (32). This 
may be helpful in addressing puzzles like the hierarchy 
problem. Secondly, one may consider non-Gaussian path 
integrals for self-interacting fermions. Finally, it would 
be interesting to study scalars coupled to gauge fields 
and, especially, to analyze the mechanism of spontaneous 
symmetry breaking using non-Gaussian path integrals. 
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